SINP-TNP/02-20 



Novel Quantum States of the Rational Calogero Models 
Without the Confining Interaction 



O ' B. Basu-Mallick 1 , Pijush K. Ghosh 2 and Kumar S. Gupta 3 

O 

^ \ Theory Group 



X 



Saha Institute of Nuclear Physics 
1/AF Bidhannagar, Calcutta - 700064, India 



Abstract 



CO 
> 

o 

We show that the A-particle A/y-i and Bn rational Calogero models without the 
harmonic interaction admit a new class of bound and scattering states. These states owe 
£S) ' their existence to the self-adjoint extensions of the corresponding Hamiltonians, labelled 

by e tz where z € R (mod 2ir). It is shown that the new states appear for all values of 
N and for specific ranges of the coupling constants. Moreover, they are shown to exist 
even in the excited sectors of the Calogero models. The self-adjoint extension generically 
breaks the classical scaling symmetry, leading to quantum mechanical scaling anomaly. 
The scaling symmetry can however be restored for certain values of the parameter z. We 
also generalize these results for many particle systems with classically scale invariant long 
range interactions in arbitrary dimensions. 
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1. Introduction 



The spectrum of the iV-particle one dimensional Aw-i rational Calogero model without 
the confining term has traditionally been described in terms of scattering states alone [1]. 
It has however recently been shown [2] that the spectrum of such a system admits negative 
energy bound states for certain values of the system parameters. These bound states owe their 
existence to the self-adjoint extension [3, 4] of the "radial" part of the Hamiltonian. These 
self-adjoint extensions are labelled by e lz where z E R. The parameter z (mod 2ir) classifies all 
possible boundary conditions for the system. Indeed, the corresponding result for the two-body 
case has been known for a long time [5], which has recently found applications in the context 
of black holes [6, 7, 8]. 

In the analysis presented in Ref. [2], the new bound states of the A^-i rational Calogero 
model were found to exist only for N = 3 and 4 and for the zero "angular momentum" 
sector of the problem. This happened due to the fact that the wavefunction was required to 
vanish when any two particles coincide. Such a boundary condition allows to smoothly extend 
the wavefunction for a given particle ordering to the whole of the N-particle configuration 
space corresponding to all possible orderings. In this paper we observe that this boundary 
condition on the wavefunction can be relaxed while still keeping it square integrable for any 
given particle ordering. The wavefunction obtained in this way can then be extended to all 
other particle orderings by using the permutation symmetry, although not in a smooth fashion. 
As a consequence of the more general boundary condition, we are now able to find the new 
bound states for arbitrary values of particle number N. Even more remarkably, the self-adjoint 
extensions are now found to exist even in the higher "angular momentum" sectors of the 
problem. 

The scattering state solutions of the A^-i model in absence of the confining term was 
already found by Calogero in Ref. [1]. The range of the coupling there was such that the 
corresponding Hamiltonian Ha n _ 1 was essentially self-adjoint. In this paper we have been able 
to extend that analysis for the parameter range where -f^A JV _ 1 is not essentially self-adjoint but 
admits self-adjoint extensions. In particular, we have found a new class of scattering states for 
each choice of the self-adjoint parameter z. The phase shifts for these scattering states depend 
explicitly on choice of z. Moreover, for generic values of z, the phase shift is found to depend 
on the momentum as well. 

The Calogero model without the confining term belongs to a class of iV-particle 

systems with long range interactions which are classically scale invariant [9-15]. The scale 
invariance manifests itself through the absence of dimensionfull coupling constants. It is thus 
natural to ask if an analysis similar to that in Ref. [2] applies to this wider class of systems, 
including the higher dimensional cases. In this paper we address this issue and show that these 
systems indeed admit a new class of bound and scattering states. These states appear due 
to the self-adjoint extensions of the corresponding Hamiltonians. As examples of this general 
formulation, we have analyzed the bound and scattering state sectors of the Aw-i and Bn 
Calogero models in one dimension. The and Cjv Calogero models in one dimension arise as 
special cases of the B N model. In all these cases we have found that the new quantum states 
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exist for arbitrary values of N and for suitable ranges of the coupling constants. Moreover, 
these states appear in the excited sectors as well. The Calogero-Marchioro model [12] in D 
dimensions has also been analyzed to illustrate our formalism in dimensions higher than one. 
It has been shown that the new states exist in this model as well with suitable restrictions on 
the system parameters. 

The classically scale invariant systems considered here are such that 0(2, 1) appears as their 
spectrum generating algebra [16, 17]. It is however known that the classical scaling symmetry 
may be broken in the presence of the self-adjoint extensions [18, 19]. We shall show that for 
the general class of systems considered here, the self-adjoint extensions typically break the 
scale invariance at the quantum level. This is evident by the appearance of the momentum 
dependence of phase shifts in the scattering sector and by the appearance of bound states 
in certain cases. The underlying reason behind this is that the domain of self-adjointness 
of the associated Hamiltonian is not kept invariant by the scaling operator [19, 20, 21, 22]. 
Scale invariance at the quantum level is recovered only for special choices of the parameter z. 
The systems considered here thus provide further examples where scaling anomaly appears in 
quantum mechanics. 

The organization of the paper is as follows. In Section 2 we provide a general analysis of 
iV-particle scale invariant systems in one and higher dimensions. Using a suitable separation 
of variables, the effective Hamiltonian in the "radial" coordinate is obtained. The quantization 
of the effective Hamiltonian together with its self-adjoint extension is discussed in Section 3. 
It is shown that new bound and scattering states emerge naturally in this analysis. In Section 
4, we discuss the application of this quantization method to several examples. These include 
An-i, Bn and Dn Calogero models in one dimension and Calogero-Marchioro model in higher 
dimensions. The restriction imposed on the system parameters necessary for the existence of 
the new quantum states is discussed in each case. We conclude the paper in Section 5 with 
some discussions. 



2. Effective Hamiltonian for Scale Invariant Potentials 



In this Section we show that the Hamiltonian of a class of many-particle systems with long 
range and classically scale invariant interactions can be related to the Hamiltonian of a single 
particle with inverse-square interaction on the half line. Let us consider a system of N identical 
particles of mass m in one dimension whose Hamiltonian in units of ImhT 2 = 1 is given by 



N 



H = T. + {d lW f - %w] , (2.1) 
1=1 

where di denotes the derivative with respect to the position of the i th particle. The function 
w appearing in Eqn. (2.1) depends on the particle coordinates and is assumed to satisfy the 
relation 

w = -lnG, (2.2) 
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where G is a homogeneous function of degree d, i.e. 



N 



J2x t d t G = dG. 



(2.3) 



i=l 



Using Eqns. (2.1) - (2.3), we can show that the many-body interaction in H scales in the same 
way as the kinetic energy term. This implies that there are no dimensionfull coupling constants 
in the system, indicating that H is classically scale invariant. It may be noted that 0(2, 1) 
appears as the spectrum generating algebra for H [16, 17, 15]. 

We would like to solve the Schrodinger's equation 



Hip = Eip. 

For this purpose, consider an ansatz of the form 

i> = e- w P k {x)ct>{r), 



(2.4) 



(2.5) 



where x = (xi, X2, ■ ■ ■ , %n) and 0(r) is taken to be a function of the "radial" variable r defined 
by r 2 = Eili xj. The function P k (x) is a homogeneous polynomial of the particle coordinates 
with degree k > 0, i.e. 



N 



satisfying the equation 



^XidiP k = kP k , 



J2d!P k + 2j2 d i wd i P k = °- 

1=1 i=l 



(2.6) 



(2.7) 



Note that for k = 0, Pq as an arbitrary constant ( which can be chosen to be unity without 
any loss of generality ) is a valid solution of Eqn. (2.7). The solutions of Eqn. (2.7) for higher 
values of k are known only for a few one dimensional many-body systems. If in a given model 
the solutions of Eqn. (2.7) for nonzero values of k are not known, we shall restrict our analysis 
only to the k = sector, with P = 1. Substituting Eqn. (2.5) in (2.4) and using Eqns. (2.3), 
(2.6) and (2.7), we get 

d 2 „ ld_ 

r dr 



dr 



0(r) = £0(r), 



where 



N 7 J 

v = 1 + k + d. 

2 



(2.8) 
(2.9) 



Lets us now define the function x(r) through the transformation 



0(r) = r 

Substituting Eqn. (2.10) in Eqn. (2.8) we get 

d 2 



-(*+") 



xir). 



H X (r) = 



dr 2 



+ 



X(r) = Ex(r). 



(2.10) 



(2.11) 



4 



The operator H in Eqn. (2.11) defines the effective Hamiltonian of the system in the "radial" 
variable. 

In deriving Eqn. (2.11) we did not make any assumption about translation invariance of 
the system. If the system is translationally invariant, then both Pk{x) and w(x) obey the same 
symmetry. In that case, we can still derive Eqn. (2.11) after eliminating the center of mass 
coordinate with r and v now being given by 

1 N 

2 1 I \ 2 



r = 



2iV . . 



Eixi-Xj) 2 , (2.12) 



TV — 1 

v = — 1 + k + d. (2.13) 

Comparing Eqns. (2.9) and (2.13), we find that these two expressions differ by a contribution 
coming from the center of mass of the system. 

In this paper we shall discuss both the scattering and bound states sectors of the Hamiltonian 
H, or equivalently, of the effective Hamiltonian H. For the purpose of discussion of the bound 
states, we now find the measure for which the eigenfunctions x( r ) would be square integrable. 
Using Eqns. (2.2), (2.5) and (2.10), we get 

V> = GP k {x)r- { ^ +u) x(r). (2.14) 

Note that r takes values in the positive real axis. For translationally invariant systems, the 
Hamiltonian H can equivalently be described using (N — 1) dimensional hyperspherical co- 
ordinates with the radius r and the (N — 2) angular variables In this case, both G and 
Pk(x) are translationally invariant homogeneous function of degree d and k respectively. The 
wavefunction ip(x) can thus be factorized as ip(x) = £(r)F(fi) where 

ar)=r d+k - v -h(r), (2.15) 

and F(Q) is a function of the angular variables f2j alone. Using the Eqns. (2.13) and (2.15), 
we see that 

* C(r)£(r)r {N - 2) dr = I ' x*(r)x(r)dr, (2.16) 

Jo 

which implies that x{ r ) £ L 2 [R + ,dr]. For translationally non-invariant systems, the Hamilto- 
nian H can be described using N dimensional hyperspherical coordinates. The same line of 
argument as above using Eqn. (2.9) instead of Eqn. (2.13) again gives x( r ) £ L 2 [R + ,dr}. 

We have so far restricted our attention where the space dimension is one. It can easily be 
shown that the above derivation is also valid in D space dimensions. In this case the particle 
coordinates would be D dimensional. Thus substituting Xj — > r { and <9j — > Vj in Eqn. (2.1) 
and following the subsequent analysis, we get the same effective Hamiltonian H with v being 
given by 

v= (N-l)D_ 1 + k + d (21?) 
for translationally invariant systems and 

ND , . 

v=— --1 + k + d (2.18) 
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for translationally non-invariant systems. Note that Eqns. (2.17) and (2.18) differ by a contri- 
bution coming from the center of mass of the system. 



3. Quantization of the Effective Hamiltonian 



In this Section we shall discuss the quantization of the effective Hamiltonian in Eqn. (2.11) 
[5]. We shall present the analysis in the scattering sector directly using von Neumann's theory 
of deficiency indices [3] when z/ 2 > 0. This will lead to an alternate derivation of the bound 
state solutions for this system as well. The operator H belongs to a general class of objects 
known as unbounded linear differential operators on a Hilbert space [3]. We shall first outline 
some properties of these operators which are needed for our purpose. 

Let T be an unbounded differential operator acting on a Hilbert space H and let D(T) be 
the domain of T. The inner product of two element a, (3 G H is denoted by (a, (3). Let D(T*) 
be the set of G H for which there is a unique i] G 7i with (T£,0) = (£,77) V £ G D(T). 
For each such G D(T*), we define T*0 = 77. T* then defines the adjoint of the operator T 
and D(T*) is the corresponding domain of the adjoint. The operator T is called symmetric or 
Hermitian iff (T<f>,r)) = ((f), Trj) V 0, r\ G D(T). The operator T is called self-adjoint iff T = T* 
andD(T) = D(T*). 

We now state the criterion to determine if a symmetric operator T is self-adjoint. For this 
purpose let us define the deficiency subspaces K± = Ker(i =f T*) and the deficiency indices 
n±(T) = dim[lf±]. Then T falls in one of the following categories: 

1) T is (essentially) self-adjoint iff = (0,0). 

2) T has self-adjoint extensions iff n + — n_. There is a one-to-one correspondence between 
self- adjoint extensions of T and unitary maps from K + into K_. 

3) If n + 7^ n_, then T has no self- adjoint extensions. 

We now return to the discussion of the effective Hamiltonian H . This is an unbounded 
differential operator defined in R + . H is a symmetric operator on the domain D(H) = {0(0) = 
0'(O) = 0, 0, 0' absolutely continuous, G L 2 (<ix)}. We would next like to determine if H is 
self-adjoint, for which we treat the cases v ^ and v — separately. 

3.1 ^ / 

The deficiency indices n± are determined by the number of square-integrable solutions of 
the equations 



respectively, where H* is the adjoint of H. Note that H* is given by the same differential 
operator as H. From dimensional considerations we see that the r.h.s. of Eqn. (3.1) should 
be multiplied with a constant with dimension of length -2 . We shall henceforth choose the 
magnitude of this constant to be unity by appropriate choice of units. The solutions of Eqn. 
(3.1) are given by 



H*<j> ± = ±i0±, 



(3.1) 



0+(r) = r^«(re^), 
0_(r) = r^H^\re-^), 



(3.2) 
(3.3) 
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where if„'s are Hankel functions [23]. The functions (f>± are bounded as r 
they behave as 



+ (r) 



0_(r) 



smz/7r 



i 

"2 e" 



1 4 



e 4 



2^ r(l + i/) 



2"" 



i 



+ 



2^ r(i + z/) 2-" r(i-i/) 



oo. When r — > 0, 

(3.4) 
(3.5) 



We see that 0± are not square integrable when v 2 > 1. In this case has deficiency indices 
(0,0) and is (essentially) self-adjoint on the domain D(H) [4]. On the other hand, both <fi± are 
square integrable when either -1 < i/ < or < i/ < 1. We therefore see that for any value 
of v in these ranges, H has deficiency indices (1, 1). In this case, H is not self-adjoint on the 
domain D(H) but admits self-adjoint extensions. The deficiency subspaces K± in this case are 
one dimensional and are spanned by the functions (f>±. The unitary maps from K + into K_ 
are parameterized by e lz where z G R (mod 2ir). The operator H is self-adjoint in the domain 
D Z (H) = D(H) © {a(0 + (r) + e lz 4>„(r))} where a is an arbitrary complex number [3, 4]. 

We would now like to obtain the spectrum of H in the domain D Z (H). We start with the 
discussion of the scattering states given by the positive energy solutions of Hamiltonian H. For 
that we set E = q 2 , where q is a real positive parameter. The general solution of Eqn. (2.11) 
can be written as 

x (r) = r-2 [a(q)J u (qr) - b(q)J_ u (qr)} (3.6) 

where a(q) and b(q) are two as yet undetermined coefficients. J u in Eqn. (3.6) refers to the 
Bessel function of order v [23]. Note that in the limit r — > 0, 



</>+(r) + e i2 0_(r) 



and 



smz/7r 



e 1 4 



2 V 



r(i + i/) 



a(q)q u 
2 V Y(l + u) 



2~ v r(i - 1/) 

b(q)q- u 



2-^ r(i-u) 



(3.7) 



(3.8) 



If x( r ) ^ D Z (H), then the coefficients of of r^ + 2 and r " + 2 in Eqns. (3.7) and (3.8) must 
match. Comparing these coefficients we get 



«(g) _ sin(f + 3tt^) 2|/ 



&(?) 



sm(§ +7ri 



(3.9) 



Next we calculate the S'-matrix and the associated phase shift. In the limit r — > oo, the leading 
term in the asymptotic expansion of x( r ) in Eqn. (3.6) is given by 



V^q 



a(q)e 



b(q)e 



+ 



-iqr 



a(q)e 



b(q)e 



(3.10) 

By dividing the coefficient of outgoing wave (e iqr ) by that of the incoming wave (e iqr ) ) we 
obtain the the S'-matrix and phase shift 5(q) as 



S{q) 



o 2iS(q) 



a(q)e~^ +1 2^ - b(q)e^-^ 



(3.11) 
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Finally, using Eqns. (3.9) and (3.11), we get 

s() = e 2*%) = <T sin(f + 37rf )e-^)f - g sin(f + Trfje^f ^ 
g-^ sin(f + 3tt| )e^ + 5)f - sin(f + tt| )e"^-5)§ ' 

Let us now consider the bound state solutions of Eqn. (2.11). For any given value of v 
in the allowed range the S'-matrix in Eqn. (3.12) has a pole on the positive imaginary axis 
of the complex g-plane. Such a pole indicates the existence of a bound state for the effective 
Hamiltonian H. By taking q = ip as the pole for the S'-matrix in Eqn. (3.12), one can easily 
derive the corresponding bound state energy E = —p 2 as 



E = 



sin(f + 3tt|) 
sin(f +7T|) j 



(3.13) 



Thus we see that for a given value of v within the allowed range, H admits a single bound state 
with energy given by Eqn. (3.13). It may be noted that for a fixed is, the bound state exists 
only for those values of z such that the quantity ^(z+n^) m Eqn. (3.13) is positive [5]. Using 
Eqns. (3.6), (3.9) and (3.13), we see that the the corresponding bound state eigenfunction is 
given by 

x {r) = Br^H^\iJ\E\r), (3.14) 



where B is the normalization constant. It may be noted that the expression of the bound state 
wavefunction for < v < 1 has been discussed in the literature [5, 6, 2] and agrees with Eqn. 
(3.14), when restricted to the same range of v. 

We would like to remind the reader that when v 2 > 1, H is essentially self- adjoint in the 
domain D(H). In this case, the scattering state solutions are given by 

Xi r ) = a(q)r^ J u (qr) , (3.15) 

which is same as Eqn. (3.6) without the term proportional to J- V {qr). No bound states exist 
in this case. 

3.2 v = 

In this case, the solutions of the equations 

H*i)± = ±itp± (3.16) 

are given by 

V>+(r) = AH^\re^), (3.17) 
V>_(r) = r^H^\re^) (3.18) 

respectively. ip± are bounded functions as r — > oo. In order to find their behaviour for small r, 
we first note that in the limit r — > 0, 

J (r) = 1 + C(r 2 ) (3.19) 
N (r) = - [7 - ln2 + lnr] + O (r 2 \nr) (3.20) 

7T \ / 
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where 7 is Euler's constant and N is the Neumann function [23]. Using Eqns. (3.17), (3.19) 
and (3.20) we see that when r — > 0, ip+(r) behaves as 



M r ) 

In the same limit, i])-{r) behaves as 



2i i, 1 
— r 2 mr + r 2 

7T 



2i 



ip-(r) — > r 2 lnr + r 



7T 



i + ^( 7 -ln2) 

Z 71 



1 %- ln2 ) 
Z 71 



(3.21) 



(3.22) 



From Eqns. (3.21) and (3.22) we see that both ijj±(r) are square integrable functions [8]. H 
therefore has deficiency indices (1,1) and the corresponding deficiency subspaces K± are again 
1-dimensional, spanned by the functions i])±{r). As before, the operator H is not self-adjoint on 
D(H) but admits a one-parameter family of self- adjoint extensions labeled by e tz where z G R. 
It is self-adjoint in the domain D Z (H) which contains all the elements of D(H) together with 
elements of the form a(ip + (r) + e lz ip_(r)) where a is an arbitrary complex number. [3, 4]. 

The scattering states associated with positive energy solutions of Eqn.(2.11) are given by 



x (r) = r 2 [a(q)J (qr) - b(q)N (qr)} , 



(3.23) 



where E = q 2 and a(q), b(q) are two as yet undetermined coefficients. In order to find the 
ratio of these two coefficients, as before we use the fact that if H has to be self-adjoint, the 
eigenfunction x( r ) m ust belong to the domain D Z (H). In the limit when r — > 0, we have 



26(g) 



T 2 In T -\- T 2 



7T 



71 



71 



(3.24) 



We now compare the coefficients of r 2 and r 2 lnr in Eqn. (3.24) with those appearing in the 
expression of e l %ip-(r) + e~ l %ip + (r). Comparing the coefficients of r 2 lnr we find 



b{q) = -2sin-. 

Comparing the coefficients of and using Eqn. (3.25) we obtain 



(3.25) 



( \ 2b ^ 1 
a(q) — In q 



71 



Z 

cos-. 
2 



Using Eqns. (3.25) and (3.26) and assuming that z 7^ 0, we finally obtain 



a(g) 
Kq) 



71 



lng 



1 z 
-cot-. 

2 2 



(3.26) 



(3.27) 



Thus we find that for any generic value of the self-adjoint parameter z 7^ 0, the ratio of a(q) 
and b(q) depends on the momentum q. We shall comment on the z = case at the end of this 
Section. 
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For calculating the S-matrix and phase shift for the above mentioned scattering process, we 
consider the leading term in the asymptotic expansion of x( r ) °f Eqn. (3.23) at r — > oo limit. 
This gives 



X(r) 



e iqr e 



[a(q)+ib(q)} + 



1 



e~ iqr e^ [a(q) - ib(q)} . 



(3.28) 



By dividing the coefficient of outgoing wave (e tqr ) by that of the incoming wave (e tqr ) and 
using Eqn. (3.28), we obtain the the S'-matrix and phase shift as 



S(q) 



;7r | laq — \ cot § + 

l~2 _7T Z Z 

- In a 

7T * 



icotf 



^3.29) 



When z ^ 0, the S'-matrix in Eqn. (3.29) again has a single pole on the positive imaginary 
axis of the complex g-plane. Such a pole naturally indicates the existence of a bound states for 
the effective Hamiltonian H . By taking q = ifi as the pole for the S'-matrix in Eqn. (3.29), we 
obtain the corresponding bound state energy E = —ji 2 as 



E = — exp 

The corresponding eigenfunction is given by 



71 Z 

— cot— 
2 2 



x (r) = Cr^K ( ^J\E\r^j = C^r^H^ (i\J\E\r^j , 



(3.30) 



(3.31) 



where C is a constant and K is the modified Bessel function [23]. 

This concludes our discussion regarding the quantization of the effective Hamiltonian H. 
We end this Section with the following general observations. 

1) We have seen that for a given value of is, each value of the parameter z (mod 2tt) produces a 
different spectrum of H. The spectrum of H, and correspondingly the parameter space of the 
quantum theory is thus characterized by the pair (u, z). Let us assume for the moment that 
v 7^ 0. Consider now two pairs of parameters (v, z) and (z/, z') corresponding to two different 
quantum theories. It is easily seen that the spectrum of these two theories are identical if 
v 1 = —v and z' = z + 2t\v. We therefore have an equivalence relation (u, z) ~ (— z/, z + 2^u) on 
the parameter space of H. The transformation (u, z) — > (z/, z') relating two different quantum 
theories is analogous to the duality symmetry in this system. In the case when v = 0, we 
automatically have (u, z) = (z/, z') for all values of z. We can thus say that the pair [y = 0, z) 
defines the self-dual points in the parameter space. 

2) As mentioned before, the classical system that we started with possesses scaling symmetry. 
However, in the presence of the self-adjoint extensions, the system admits bound state(s) and the 
phase shifts in the scattering sector depend explicitly on the momentum. These are indicative 
of the breakdown of the scaling symmetry at the quantum level. We first analyze this issue 
when — 1 < v ^ < 1. Let us consider the action of the scaling operator A = -^(r-^: + J^r) on 
an element 0(r) = 0+(r) + e* 2 0_(r) G D Z (H). In the limit r — > 0, we have 



A(f)(r) 



smz/7r 



2^ 



r(i + i/) 



~2~» r(i - v) 



(3.32) 
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In order for A0(r) G D Z (H), we must have A0(r) ~ C<p(r) where C is a constant. However, 
the two terms on the r.h.s. of Eqn. (3.32) are multiplied by two different factors , i.e. (1 + v) 
and (1 — v). Due to the presence of these different multiplying factors, we see that A0(r) in 
general does not belong to D Z {H). Scale invariance is thus broken at the quantum level for 
generic values of z. However, from Eqn. (3.32) it is clear that for special choice of z = — i y, 
A0(r) G D Z (H) and the scaling symmetry is recovered [6]. In addition, we find that the scaling 
symmetry is also preserved at the quantum level when z = — For these choices of z, the 
bound states do not exist and the S matrix becomes independent of the momentum. The 
scaling symmetry thus is anomalously broken due to the quantization for generic values of z. 
In the case when v — 0, a similar analysis as above again shows that the self-adjoint extension 
generically breaks the scaling symmetry. For this case, the scale invariance can be recovered at 
the quantum level only for z = 0. 



4. Examples 



We have shown in the previous Section that the effective Hamiltonian H, and thus the 
Hamiltonian H, admits self-adjoint extensions when — 1 < v < 1. As discussed before, v is a 
function of several system parameters. The condition on v implies that H admits self- adjoint 
extensions only if the system parameters are suitably constrained. In this Section we shall 
consider several models with suitable ranges of parameters, such that the corresponding Hamil- 
tonians admit self-adjoint extensions. The new bound and scattering states will be obtained 
by substituting the solutions of eigenvalue Eqn. (2.11) into Eqn. (2.14). As mentioned before, 
the effective Hamiltonian is essentially self- adjoint when v 2 > 1. This is the case considered by 
Calogero in Ref. [1] for the A N _i model. In the examples discussed below we shall not consider 
the range v 2 > 1. 

4.1 A Ar _i Calogero Model 

The Hamiltonian of the Calogero Model without the harmonic term is given by 

Ha n _, = - E S-2 + (« 2 - \) E 0* - *;)" 2 > (4-1) 

where i,j = 1,2, ■■■,N. This Hamiltonian can be obtained from Eqns. (2.1) and (2.2) by 
choosing G as 

G = \[{x l -x ] y + K (4.2) 

i<j 

which is a homogeneous function of degree d = (a + ^)N(N — l)/2. The Hamiltonian Ha n _ 1 is 
translationally invariant and the parameter v is determined from Eqn. (2.13) as 

, / 1\N(N-1) N — 3 
v = k + {a + -)-L- r ->+—. (4.3) 

For this model, it is known that Eq. (2.7) determining Pk{x) is exactly solvable for any k > 
[1, 9]. Below we shall discuss both bound and scattering state sectors of the Hamiltonian H An i 
for arbitrary k > 0. 
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The bound state sector of this system was analyzed in Ref. [2] under the assumption 
that a + \ > 0. In this case, the wavefunction vanishes in the limit when any two particle 
coordinates coincide, though the corresponding current might show a divergent behaviour. 
Such a boundary condition on the wavefunction is quite similar to what one encounters in the 
case of (^-function Bose gas with infinitely large value of two-body coupling constant [10, 24]. 
This type of boundary condition allows one to construct continuous eigenfunctions on the whole 
of configuration space by first solving the eigenvalue problem for a definite ordering of particles 
e.g. x\ > x 2 > • • • > xn and then smoothly extending it to the rest of the configuration space 
using permutation symmetry associated with identical particles. It is obvious from Eqn. (4.3) 
that for a + | > and iV > 3, v is a positive definite quantity. Thus, the allowed ranges of 
v for H An i to have bound states is further restricted to < v < 1. Analyzing Eqn. (4.3) 
with the constraints on a + | and v discussed above, it was shown in Ref. [2] that the A N _i 
Calogero model admits negative energy bound states only for N = 3,4 and k = 0. 

The bound state wavefunctions of Ha n _ 1 obtained in Ref. [2] are normalizable, although 

N — 3 

as r — > 0, they have a singularity of the form r~ v ~ . Note that the singularity at r = 
corresponds to the case where all the particles coincide at the same point. This is different 
from the singularity arising from the coincidence of any two particles, which is avoided by 
imposing the constraint a + | > 0. However, from the viewpoint of self-adjoint extensions 
alone, there are no reasons a priori to distinguish between the singularity at r — > limit and 
the singularity arising from the coincidence of any two particle coordinates. In view of this, we 
now let both these types of singularities to appear in the wavefunction. It may be noted that for 
both bound and scattering state solutions, the angular part of the total wavefunction must be 
square integrable. The angular part of the wavefunction receives a contribution from the factor 
G in Eqn. (4.2), which is singular when a + \ < 0. Requiring square integrability of the angular 
part of the wavefunction puts the restriction that a+| > — \. A negative value for the parameter 
a+| in this range leads to a singularity in the wavefunction resulting from the coincidence of any 
two particle coordinates. This restricts the range of a continuous eigenfunction within a region 
of configuration space corresponding to definite ordering of particles like X\ > x 2 > • • • > x^. 
Using permutation symmetry, such an eigenfunction can be extended to the whole of phase 
space, although not in a smooth fashion. The bound state wavefunctions of Ha n1 thus obtained 
are completely normalizable for a + \ > —\. 

We now analyze the case where N is fixed while a + \ is allowed to vary subject to the 
restriction that a + \ > —\. We note that since — 1 < v < 1, we must have, 

N- 1 + 2k 1 A- 1 + 2k 4 ,„ „, 

<a + -< — — — 7^ + ^7777 7T- (4.4) 



N(N-l) 2 N(N-l) N(N-l)' 

Imposing the condition that the upper bound on a + | in Eqn. (4.4) should be greater than 
— |, we find that k is restricted as 

k < i (iV 2 - 3N + 10) . (4.5) 

Thus, the allowed values of k are 0, 1, 2, . . . , K = {(N 2 - + 10) /4}. The symbol {x} denotes 
the integral part of x if x is non-integer and is equal to x — 1 for integer x. For a fixed value 
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of N, the self-adjoint extension is admissible when the coupling constant lies within the range 



-2 <a+ 2 < N(N^Ty (46) 

Note that the upper bound in Eqn. (4.6) is negative for N > 5. This implies that the interaction 
in Eqn. (4.1) is repulsive for N > 5. 

We have seen in the previous paragraph that for fixed value of N and variable coupling, the 
self-adjoint extension and, hence, new quantum states are admissible for several values of k. We 
now ask as to what are the allowed values of k when both N and a + \ are kept fixed. To this 
end, we introduce a quantity [3{a, N) = (a + \)N(N - 1)/2 + (N - 3)/2 so that v = k + (3{a,N). 
Note that for fixed values of a and N, f3 is fixed. Moreover, since a+|>— |,/5is bounded 
from below, f3 > —f3 Q = —\(N 2 — 3N + 6). It is now easy to see that with — 1 < v < 1 and 
k > 0, number of allowed values of k for fixed (3 is at most 2. In order to find the specific values 
of k, we consider the cases v ^ and v — separately. 

4.1.1 v / 

We now use the constraint, — 1 < u ^ < 1, and present our results in the table below. 



Bounds on (5 


Ranges of a + | 


Allowed values of 


< p < 1 
-1< (3 < 
-2 < (3 < -1 

-(Z + l) <(3<-l 

-{K-l) <f3< -(K-2) 
-(3 <(3< —(K - 1) 


, 1 > 5-iV 
u ^ 2 — N(N—1) 

3-N < . 1 < 5-7V 
N(N-l) ^ 1 2 ^ N(N-l) 

^ ^ u -r 2 v. N(N—1) 
N(N-l) ^ " ^ 2 ^ AT 

AT+2/-1 x , 1 / 7V+22-3 
N(N-l) ^ 1 2 ^ N(N-l) 

N+2K-5 ^ , 1 / N+2K-7 
N(N-l) ^ 1 2 ^ N(N-l) 

1 / „ , 1 / JV+2K-5 

2 ^ u -T 2 ^ N(N-l) 


X 


0, 1 
1,2 

U + i 

K-2, K-l 
AT-1, AT 



It is curious to observe how only two successive values of fc are allowed within a specific range 
of a + |, although any lower value than this within the same range is not allowed. 

When v ^ 0, using Eqns. (2.14), (3.14) and (4.2), we see that the bound state wavefunction 
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for Ha ni is obtained as 

^ = BX\( Xl - x ) a+ ^P k (x)r-»H^\iJ\E\r). (4.7) 

i<j 

The corresponding bound state energy E is given in Eqn. (3.13). Note that for fixed values of 
N, a and z, the system may admit two bound states corresponding to two possible values of 
k. Similarly, the scattering state solutions of Ha n-1 obtained by using Eqns. (2.14), (3.6) and 
(4.2) are given by 

1> = I]>* - x^Puixy-" [a(q)J u (qr) - b(q)J_ u (qr)} , (4.8) 

i<j 

where ^j- is given in Eqn. (3.9). 
4.1.2 v = 

We now discuss the case for which v = 0, or equivalently f3 = —k. It follows from Eqn. (4.3) 
thata + | = ^pEfy. Since a + \ > — §, the allowed values of A; are 0, 1, ... , {(iV 2 — 3iV + 6)/4}. 
In this case, using Eqns. (2.14), (3.31) and (4.2) we see that the bound state wavefunction for 
the Hamiltonian Ha n _ 1 is given by 

V> = C T lite - x j )'**P k {x)r- t 'HW (i\[\E\r) , (4.9) 
z i<j v / 

where the bound state energy E are given in Eqn. (3.30). Similarly, the scattering state 
solutions in this case are obtained using Eqns. (2.14, (3.23) and (4.2) as 

^ = IIO* - Xj)^P k (x)r- v [a(q)J (qr) - b(q)N (qr)} , (4.10) 

i<j 

where is given in Eqn. (3.27). 

We end this Section with the following observations: 

(i) Self-adjoint extension and, hence, the existence of the new quantum states is admissible in 
the Am -i Calogero model for arbitrary N. 

(ii) The self-adjoint extension for k > is remarkable in the following sense. If a Hamiltonian 
admits self-adjoint extension, it usually occurs in the zero angular momentum sector of the 
Hamiltonian. We are not aware of any counter-example to this fact. For the if J 4 JV _ 1 , k can 
be related to the angular-momentum eigenvalue of the Hamiltonian [11]. Considering the case 
k 7^ amounts to studying the model in the non-zero angular momentum sector of the model. 
Thus, we provide the first example in the literature where self-adjoint extension is admissible 
in non-zero angular momentum sector of a Hamiltonian. 

4.2 B N Calogero model 

The Hamiltonian of the Calogero model without the harmonic term is given by 

N 

h b n = -Ei^ + ^i- 1 )!^-^)" 2 + (x i + x j y 2 

i=l aX i ift 
N 

+92(92 - 1)E^" 2 - (4.H) 
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The Hamiltonian in Eqn. (4.11) can be obtained from Eqns. (2.1) and (2.2) by choosing G as 



i<j s 

where s — 1, 2, • • • , N. G in Eqn. (4.12) is a homogeneous function of degree d = g\N{N — 
1) + g 2 N. It is known that for the B N Calogero model, Eqn. (2.7) admits solutions only for 
even k > [9]. Moreover, the Hamiltonian Hb n is not translationally invariant. Consequently, 
v is determined from Eqn. (2.9) as 

v = y - 1 + 9l N(N - 1) + g 2 N + 2k, (4.13) 

where k — 0, 1, 2, Using the constraint — 1 < v < 1, it can be shown that for fixed values 

of N, <7i and g 2 , there is at most one allowed value of k. 

We encounter three kinds of singularities in the Hamiltonian (4.11). They occur whenever, 

(i) any two particles coincide, 

(ii) any particle is at the position of the image of any other particle, 1 • 6 • X i ^ X >i - 

(iii) any particle is at the origin, i.e. x-i — > 0. 

If we do not want the wavefunction to reflect singularities of these kinds, we have to put the 
restriction gi,g 2 > 0. For such cases, the new states exists only in the k = sector and for 
N = 2 and 3. Compared to the similar situation in the case of Calogero model, N = 4 is 

not allowed for the present case. This may be attributed to the fact that Bn Calogero model 
is not translationally invariant. 

In the more general case, following the discussion for the A^^i Calogero model, we allow 
the wavefunction to have singularities while maintaining the square-integrability of the angular 
part, which demands that gi,g 2 > — §■ We would now like to find the allowed values of k when 
N is kept fixed but the couplings are allowed to vary. To this end, we first write g\ and g 2 as, 

gi = ~\ + el g 2 = ~ + e 2 2 , (4.14) 

so that ei and e 2 can take any nonzero real values. Using the constraint — 1 < v < 1, we find 
that 

1 2k 2 t\ 1 2k 2 



_ < e H - < I (4 15) 

2 N(N-l) 1 N-l 2 N(N-l) N(N-l)' v ' ; 

The upper bound in Eq. (4.15) restricts k as 

k <^N{N + (4.16) 

The allowed values of k are given by k = 0, 1, 2, . . . , K = {\N(N - 1) + 1}. When k = K, the 
l.h.s. of Eqn. (4.15) becomes negative if \N(N — 1) is not an integer. In that case we replace 
the lower bound in Eqn. (4.15) by 0. It is evident from the above analysis that for a fixed value 
of N, new states exist for the range of the couplings given by 

< £ ? + iv^n 4 + Niihry (417) 
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which defines an elliptical region in the €\ — e 2 plane. Since g±,g 2 > — \, the lines given by 
ei = and e 2 = are excluded from this elliptic region. For different allowed values of k, this 
region naturally separates into disjoint elliptical shells defined by Eqn. (4.15). The outermost 
shell corresponds to k = while the innermost one corresponds to k — K. There are no new 
quantum states for the values of the couplings corresponding to the boundary separating two 
consecutive shells. 

When v — 0, the couplings satisfy the relation 

e^ + ^L = -+ !~ 2A \ . (4.18) 
1 N-l 2 N(N-l) V ; 

The allowed values of k in this case are found to be k — 0, 1, ... , {(N 2 — N + 2)/4}. In general, 
a shell defined by Eqn. (4.15) contains an ellipse of the form (4.18) with the same value of 
k. However, when \N(N — 1) is an integer, the value of v — cannot be achieved within the 
innermost shell corresponding to k — K. 

When v 7^ 0, the bound state wavefunction of the B N model is given by 



i> = B\{{x 2 -x^Wxf P 2k {x)r- u H£\i^j\E\r), (4.19) 

i<j s 

where the bound state energy E is given by Eqn. (3.13). P appearing in Eqn. (4.19) is a 
solution of Eqn. (2.7) for the B N Calogero model. Similarly, the scattering state solutions are 
given by 

1> = Ife 2 - x )) 91 U x s 2p ™( x >~" - Kq)J- v (qr)] , (4.20) 

i<j s 

where |^ is given by Eqn. (3.9). When v — 0, the bound state wavefunction is given by 

^ = Cy \[{x 2 - x 2 r\[xfP 2k {x)r- v H^ (iy/\E\r} , (4.21) 

where the bound state energy E are given in Eqn. (3.30). Similarly, the scattering state 
solutions are given by 

ij = ]J(x 2 - x)Y^xfP 2k {x)r- v [a(q)J (qr) - b(q)N (qr)} , (4.22) 

i<j s 

where |||y is given by Eqn. (3.27). 

4.2.1 D N Calogero Model 

The Dn Calogero model can be obtained from the model by putting g 2 = 0, i.e. e 2 = \. 
Using the expression for v in (4.13) and the constraint — 1 < v < 1, we find, 

N + Ak N + Ak 2 . Ano . 

<9l< — 77T77T7 TT + TTTTr 7T- ( 4 - 2 3) 



2N(N-1) 2N(N-1) N(N-l) 

The constraint g\> —\ restricts k as 

1 



k < 1(N 2 -2N + A). (4.24) 
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Thus, the quantum number k can take any integral values, k — 0, 1, 2, . . . , K = {\{N 2 — 2N + 
4)}. It is clear that for a fixed value of N, new states exist for the range of g\ given by 

1 4 - N 

-2 <3l< 2N(N-iy (425) 

Note, however, that the points given by g 1 = —jw^T) belonging to the range in Eqn. (4.25) 
are excluded. For fixed values of g 1 and N, the number of allowed values of k is at most one. 
When v ^ 0, the bound and scattering state wavefunction in this case are obtained from Eqns. 
(4.19) and (4.20) respectively by putting g 2 = in them. When v — 0, g\ — — 2n{n-i) wr th 
k — 0,1, ... , {{N 2 — N + 2)/4}. The bound and scattering states wavefunctions in this case are 
obtained from Eqns. (4.21) and (4.22) respectively by putting g<i = in them. 

Finally, as in the case of A^-i Calogero model, we have the remarkable results that for 
suitable ranges of couplings, the B^ and Dn Calogero models admit self- adjoint extensions 
and consequently, new scattering and bound states for arbitrary N . These new states exist 
for certain non-zero angular momentum sector also. For fixed values of iV and the relevant 
couplings, i^Ajv-i admits a maximum of two bound states corresponding to at most two possible 
values of k. On the other hand, under similar conditions, and Calogero models admit 
only one bound state. It may be noted that the Hamiltonian of the Cn Calogero model is 
identical to Hb n except that the last term in Eqn. (4.11)is replaced by ^3(^3 — 1) X^i(2xj) -2 
where g 3 is a constant. The analysis presented above can be extended to the Cn Calogero 
model as well. 

4.3 Calogero-Marchioro Model : A D Dimensional Example 

All the models that have been discussed so far are one dimensional. We now give an 
example of a D dimensional many-particle system, known as the Calogero-Marchioro model in 
the literature [12, 13]. The Hamiltonian for this model is given by 

N 1 

ffcM = -E^ + 3(9 + 5-2)ET+9 2 E -&> (4-26) 

i=l v£j ij i¥=j¥=k V lk 

where = — r^. The Hamiltonian in Eqn. (4.26) can be obtained from D-dimensional 
analogue of Eqn. (2.1) by choosing 

G = U4 ( 4 - 27 ) 

G is a homogeneous function of degree d = gN(N — l)/2. H C m contains both two-body and 
three-body interaction terms. For D — 1, the three-body interaction term vanishes identically 
and H C m reduces to the the H An _ 1 . For D > 2, although infinitely many exact solutions corre- 
sponding to the radial excitations are known, the model is not exactly solvable. In particular, 
Eq. (2.7) determining solutions for is not known for k > 0. We will henceforth restrict 
our attention to the k = sector of the model. It may be mentioned here that the Calogero- 
Marchioro model in D = 2 is related to several systems of physical interest, e.g. quantum 
Hall effect, quantum dots, random matrix theory etc. [13]. Moreover, Hcm in D = 2 can be 
embedded into a Hamiltonian with extended Af = 2 superconformal symmetry [13]. 
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The Hamiltonian in Eqn. (4.26) is translationally invariant and the parameter v is deter- 
mined from Eqn. (2.17) as 

u = -{N -1)D -1 + -N{N -1). (4.28) 

The many-body interactions are singular at the coincident positions of any two particles. If we 
do not want the wavefunction to have singularities of this kind, we have to impose the condition 
g > 0, which requires that < v < 1. Under this condition, for fixed values of N and D, Hqm 
admits self-adjoint extensions when 

1 + ^<JV<1 + 1 (4.29) 

First note that the result of the An-i Calogero model is reproduced for D = 1. For D = 2, 
only N = 2 is allowed and for D > 3, there are no valid solution. 

In the more general case, we allow the wavefunctions to have singularities while maintaining 
the square-integrability of the angular part coming from G. This puts the restriction that 
g > — \. In this case v satisfies the constraint —\<v< \ which restricts g as 

D DA , 

-n < 9 < -n + nw^ty (430) 

Let us first consider the case v ^ 0. The condition g > — \ is satisfied for 

N 4 

D<— + — — . (4.31) 
2 N- 1 y ! 

The self-adjoint extension and consequently, the existence of scattering and bound states are 
admissible, provided Eqns. (4.30) and (4.31) are satisfied simultaneously. For D > 3, there 
exists a lower bound on N in order to have the self-adjoint extension. However, interestingly 
enough, for the physically important cases of D = 1,2 and 3, the self-adjoint extensions are 
admissible for any arbitrary N. 

When v — 0, the coupling g takes the value given by 
The constraint g > — | leads to the condition 

D< ih + T < 4 - 33 > 

Note that Eqn. (4.33) is satisfied for arbitrary N when D = 1,2. For D — 3, it is valid for 
iV > 6. We end this discussion by noting that the wavefunctions and the spectrum for this 
model can be written explicitly as before, which we do not discuss here in detail. 
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5. Conclusion 



In this paper we have discussed the quantization of N particle systems with classically scale 
invariant long range interactions. The effective Hamiltonian for these systems in the "radial" 
direction contains an inverse square term whose coefficient {u 2 — |) depends on the particle 
number N, the couplings and the "generalized angular momentum" k. It has been found that 
when — 1 < v < 1, the Hamiltonians for this entire class of systems admit a one parameter 
family of self-adjoint extensions labeled by e lz where z G R (mod 2n). The parameter z 
classifies all possible boundary conditions for which the Hamiltonian is self-adjoint. Moreover, 
the spectrum of the Hamiltonian depends explicitly on z. Each choice of the parameter z thus 
gives rise to an inequivalent quantization of the system. 

We have illustrated the general approach through several examples which have been studied 
in detail. The examples in one dimension include and Calogero models without the 

confining term. The D N Calogero model has been studied as a special case of the B N case. In 
these cases we have found that a new class of bound and scattering states appear for arbitrary 
values of the particle number N and within certain ranges of the coupling constants. Moreover, 
these new states appear not only in the k — sector, but for higher values of the quantum 
number k as well. To our knowledge, this is the first demonstration of the existence of self- 
adjoint extensions in the excited sectors of a system. This result may be attributed to the 
highly correlated nature of the many-body interaction. It may be mentioned that there is an 
important difference between the spectrum of H An x and H Bn . For fixed values of N and the 
relevant couplings, Ha n _ 1 admits a maximum of two bound states corresponding to two allowed 
values of k. On the other hand, under similar conditions, Bn Calogero models admit only one 
bound state. This can be attributed to the fact that i^A^-i is translationally invariant, while 
H Bn is not. 

As an example of our general discussion in dimension higher than one, we have analyzed the 
Calogero-Marchioro model in D dimensions. In this case, the solution of the angular equations 
is known only for k — and we have presented the analysis only for this sector. For the 
physically important cases of D — 1,2 and 3, the new states are found to exist for arbitrary 
number of particles N and within certain range of the coupling constant. 

Although the class of systems considered here is classically scale invariant, we have found 
that in presence of the self-adjoint extension, the systems may admit bound states. Moreover, 
the associated S matrix and the phase shifts are found to depend explicitly on the momen- 
tum. This happens as the classical scale invariance is broken due to quantization, which is 
manifest by the fact that the scaling operator does not leave the domain of self-adjointness 
of the Hamiltonian invariant. The analysis here thus provides further examples of quantum 
mechanical scaling anomaly. We have also shown that scale invariance at the quantum level 
can be implemented for special choices of the parameter z. It may also be mentioned in this 
context that the effective Hamiltonian H in Eqn. (2.11) can be shown to be related to the 
Virasoro algebra [25] with central charge c = 1 [2, 7]. The systems studied here would also be 
related to the Virasoro algebra. 

We have restricted our attention to the case where the coefficient of the inverse square 
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potential is not too strongly negative in order to avoid the "fall to the center" [26]. The 
strong coupling region of the analogous 2-body problem has been analyzed in the literature 
using renormalization group techniques [27]. It would be interesting to consider the problem of 
renormalization in the presence of the self-adjoint extension. The N-particle rational Calogero 
models are also related to black holes [28] and Yang-Mills theories [10, 29]. It is plausible that 
the quantum states found here would have analogues in those well. It would also be 

interesting to investigate the Calogero models associated with the exceptional Lie groups [9]. 
Finally we would like to mention that an analysis similar to the one presented in this paper 
can be done in the case of rational Calogero models in the presence of the confining term as 
well [30]. 
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